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CHAPTER 1 


The Hoyle-Narliker ‘theory 


of Gravitation 


4. introduction 


The success of Maxwell's equations has led to 
electrodynamics being normally formulated in terms of fields 
that have degrees of freedom independent of the particles in 
them. However, Gauss suggested that an action-at-a-distance 
tneory in which the action travelled at a finite velocity 
might be possible. ‘his idea was developed by wheeler and 
Feynman (1,2) who derived their theory from en action-principle 
that involved only direct interactions between pairs of part- 
icles. A feature of this theory was that the 'pseudo'-fields 
introduced are the half-retarded plus half-advenced fields 
claculated from the world-lines of the particles. However, 


2N 
Wheeler and Feynman, and, in a different way, Hogarth (3) 


were able to show that, provided certain cosmological 


コー 


conditions were satisfied, these fields could combine to 
A ; " i 4 T 
give the observed field. Hoyle and Narlikar (4) extended tne 
tueory to general space-times and obtained similar theories 
eo | qa "4 (5) } TT i {e i と I (6) 
for their 'C'-field and for the gravitational field è 


It is with these theories that this chapter is concerned. 


It will be shown that in an exp.ndinz universe tne 


advanced fields are infinite,and the retarded fields finite. 
Tais is because, unlike electric charges, all masses have tho 


same sign. 


2. The Boundary Condition 
Hoyle and Narlikar derive their theory from the 


action: 


A - £Z((G(a,b)dadb, 


Q +b 
where the integration is over the world-lines of oarticles 
a, ME In tbis expression G is a Green function 
that satisfies the wave equation: 
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where g is the determinant of Qe . oince the double sum 
in tne action A is symmetrical between all pairs of 
particles a,b , only that part of G (a,b) that is 


symmetrical between へ and b will contribute to the action 


i.e. the action can be written 


A è» Ez INCOLE: b 
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Thus ( must be the time-symmetric Green function, and can 


where G^ Co b) - 


Ko nd L 
be written: f = 2 Get 13 adv where Get 


functions. 


and lody she retarded and advanced Green 
By requiring that the action be stationary under variations 
of the RET y Hoyle and Narlikar obtain the fielá-equations 
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where mx) = ine “(x a) da . However, as a 


the particular choice of Green function, 


consequence of the 
contraction of the field-equations is satisfied identicall 


‘here ave thus only 9 equations for the 10 components of Gi; 
end the system is indeterminate. 


M =const., 


Narlikar therefore impose fm = o 


ioyle and 
as the tenth equation. By then makin. the 'smooth-fluid' 
| 5 (a) (b) 
by putting bi a mm Lie A 
+b 
they obtain the Einstein field-equations: 
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approximation, that is 
J 


Thereis an important difference, however, between these 


field-equations in the direct-particle interaction treory 
the usual general theory of relativity. In the 
metric that satisfies the 


and in 
general theory of relativity, any 


f. 7 


the field-equations is admissible, but in the direct-particle 


interaction theory only も nose solutions of the field-equation 


are admissible that satisfy the additional requirement: 


mol). £n Mn): £(G 6s, 0) da 


ilf Gurl ada + 1 £ (0,5, (2,0)da 


This requirement is highly restrictive; it will be shown 
thst it is not satisfied for the cosmological solutions of 
the Kinstein field-equations, and it appears thst it cannot 
be satisfied for any models of the univer.e thot either 
contain an infinite amount of matter or undergo infinite 
expansion, 
The difficulty is similar to that occurring in 

Newtonian theory when it is recognized that the universe 


might be infinite. 
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The Newtonian potential Q obeys the 


od - - KP (p »0) 


where È is the density. 
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Tn an infinite static universe, P would 


Pn 
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the source always has the same sign. 
ved when it was realized th:t the universe 
in 


equation is finite by a sort of'red-shift' 


be infinite, since 


difficulty 


was expending, since 


an expanding universe the retarded solution of the above 


effect. ‘The 


advanced solution will be infinite by a 'blue-shift' 


effect. 


This is unimportant in Newtonian theory, since one is free 


to choose the solution of the equation and so mey ignore the 


infinite advanced solution and take simply the finite 


retarded solution. 


similarly in the direct-particle interaction theory the 


Mm -field satisfies the equation: 


Om + ¿Rm soU (N20) 


where AJ is the density of world-lines of 


particles. 


AS in 


the Newtonian case, one may expect that the effect of the 


expansion of the universe will be to make 


the retarded solution 


finite and the advanced solution infinite. However, 


now not free to choose the finite retarded solution, 


one is 


for the 


equation is derived from a úirect-particle interaction action- 


principle symmetric between pairs of particles, and one must 


choose for M half the sum of the retarded and advanced 


solutions. We would expect tais to be infinite, and this is 


Shown to be so in the next section. 


5. The Cosmological solutions 


he Robertson-Walker cosmological metrics have the 


form yo TA n yy " 20d $3)) 
as ud a R" (E) [dr y $ db + Sin ad | 


Kr ^ 


Since tney are conformal!y flat, one can choose coordinates 


in which they become 


ds? " d Fa? - dp + pad * P 'eaigd の 2) 
a, des e x 
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wüere Job is the flat-space metric tensor and 
^. 
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For example, for the Hinstein-de Sitter universe 
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For the steady-state (de Sitter) universe 
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The Green function Y (a, b) obeys the 


equation 
GG, b) + ERG Gb) = S lab) 
ー 0 
E Es it follows tnat 


E i ab 9 "ab ; e 
d Ae QU gs He (shi a^ 


If we let n: AUR 


Ad (rias) Pn 
X 


imply the flat-space Green function equation, 


anci. 


nence ; "T ca N 
TEME Af.) に è.) 
Er LE SL(t)f 
+ glp + n) 
The'm -field is given by | PE i, ei a 
(る ) 2b "ルソー dx" oe $ A3 TM 


For universes without creation (e 
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g. the Einstein-de Sitt 


rv = 


const. For 


universes with creation (steady state) IN. fL fiz .5 const. 
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where the integration is over the future light cone. This 
will normally be infinite in an expanding universe, e.g. in 
the ijinstein-de sitter universe. 
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in the steady-state universe 
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By contrast, on the other hand, we have 


NO yrr dt 
Met (a) * A Di) 


where the integration is over the past light cone. This will 


normally be finite, e.g. in the Linstein-de sitter univers: 
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while in the steady-state universe 
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Thus it can be seen that the solution M = const. of the 


equation 


l 
wai TE 


is not, in a cosmological metric, the half-advonced plus 
half-retarded solution since this would be infinite. In fact 
in the case of the Einstein-de bitter and steady-state metrics, 


it is the pure retarded solution. 


+. The 'C'-iielà 


… — 


Hoyle and Narlikar derive their direct-particle 
/ 


interaction theory of the 'C'-field from the action 


Á 。 £ £ (f Gla, Due da db" 
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where the suffixes t, b refer to differentiation of 


^ 
G (a,b) on the world-lines of G, b resvectively. 
AN 


G is a Green function obeying the equation 


Erg CX, a 2 á'OG X). 
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we define the 'C'-field by ‘ 
レ 


Cia) a ACTION da , 
and the matter-current F by | 
14) = 2 9 Cy.b)A5 | 
Then C (x) 4 rale, 3 JI), " /-9 d : 
LIC X! jK 


ve thus see that the sources of the 'C'-fiela are the places 


where matter is created or destroyed. 


As in the case of the 'M'-field, the Green function 


must be time-symmetric, that is A ) 


C (a,b) s E Cpa, hdr TAA b 


Hoyle and Narlikar claim thet if the action of the 


1C'-fie1G is included alon;; with the action of the '™ '-field, 
a universe will be obtained that approximates to the steady- 
state universe on a large scale although there may be local 
irregularities. In this universe, the value of C will be 
finite and its gradient time-like and of unit magnitude. 

Given this universe, we may check it for consistency 0y 
claculating the advanced and retarded 'C'-fields and finding 
if their sum is finite. we shall not do this directly but 
will show that the advanced field is infinite while the retar- 
ded field is finite. 


Consider a region in space-time bounded by & three- 


dimensional space-like hypersurface D at the present time, 
and tne past light cone 5 of some point P to the future 


of Y. 


By Gauss's theorem 


ロ C/ーq dx" = dS 


SCUTIS 


Let the advanced field produced by sources within YA be L > 


nl / 、 
Then C and d C will be zero on PA 。 and hence 
dn 


Th V-gsdx* « f dS. 


s K 
But d | is the rate of creation of matter= AÀ (const.) in 
i AK 


the steady-state universe, and hence 


g S < nV. 
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p^ 
の Jn 


bl A " DAR 3 
as the point F is taken further into the future, the volume 
of the region V/ tends to infinity. However, the area of the 


A 
hypersurface D tends to à finite limit owing to horizon 
の LI 
effects. herefore the gradient peso must be infinite. 
の ん 
A similar calculation shows the gradient of the retarded 


field to be finite. ‘their sums cannot therefore ive the 

field of unit gradient required by the Hoyle-Narlikar theorr. 
It is worth noting that this result was obtained 

without assumptions of a smooth distribution of matter or 


な 


Ol 


conformal.iflatness. 


5. Conclusion 


It is one of the weaknesses of the Einstein theory of 
relativity that although it furnishes field equations it does 
not provide boundary conditions for them. ‘hus it does not 
give a unique model for the universe but allows a whole series 
of models. Clearly a theory thet provided boundary conaitions 
and tous restricted the possible solutions would be very 
attractive. The Hoyle-Narlikar theory does just thst(tne 
requirement thet NM = if. ol * i^v is 

equivalent to a boundary condition). Unfortunately, as 

have seen above, this condition excludes those models that 


seem CO correspond to the actual universe, namely the 


Robertson -wallker models. 


The calculations given above have considered the universe 
as being filled wita a uniform distribution of matter. "nis 


is legitimate if we are able to make the 'smooth-fluid' 
approximation to obtain the Hinstein equations. Alternatively 
if this approximation is invalid, it cannot be said that the 
theory yields the Einstein equations. 
It might possibly be that local irregularities could make 
Mady finite, but this has certainly not been demonstrated 
and seems unlikely in view of the fact that, in the Hoyle- 


C'-fielà, 


Narlikar direct-particle interaction theory of their 


which is derived from a very similar action-principle, it can 


be shown without assuming a smooth distribution that the 
advanced 'C' field will be infinite in an expanding universe 
with creation. 

The reason that it is oossible to formulate a direct- 
particle interaction theory of electrodynamics that does not 
encounter this difficulty of having the advanced solution 
infinite is that in electrodynamics there are equal numbers 
of sources of positive and negative sign. Their fields can 
cancel each other out and the total field can be zero apart 
from local irregularities. This sugzest that & possible way 
to save the Hoyle-Narlikar theory would be to allow messes of 
both positive and negative sign. The action would be 

e * Y d dy 
Aa EE 9。 (he (a, b) da clb (qus =! 


ad 
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where Ta, LL are gravitational charges analo;:ous to 


electric char:es. Particles of positive G in a positive 


'm '-field and particles of negative q. in a negative 'm '- 
field would have the normal gravitational properties, that is, 


toey would have positive gravitational and inertial masses, 


A particle of negative qQ in a positive 'm'-fielà would 


~ 


still follow a geodesic. Therefore it would be attracted 
a particle of positive a. . Its own gravitational effect 


however would be to repel all other particles. ‘thus it would 


^7 
Ce 
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have the properties of the negative mass described by Bondi 
that is, negative gravitational mass and negative inertial 
mass. 

Since there does not seem to be any matter having 
tnese oroperties in our region of space ( wherem const. > QO) 
there nust clearly be separation on a very large scale. It 
would not be possible to identify particles of negative の 
with antimatter, since it is known that antimatter has positive 


inertial mass. Uwever, the introduction of negative masses 


would probably raise more difficulties than it wuld solve. 
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PERTURBATIONS 


1. Introduction 
Perturbations of a spatially isotropic and homogeneous expanding 
universe have been investigated in a Newtonian approximation by 
E oU) ans relativistically by Lifshitz ^), Lirtehitz ana 
E atnikov 7 and Irvine), Their method was to consider small 
variations of the metric tensor. This has the disadvantage that the 
metric tensor is not a physically significant quantity, since one 
cannot directly measure it, but only its second derivatives. It is 
thus not obvious what the physical interpretation of a given 
perturbation of the metric is. Indeed it need have no physical 
significance at all, but merely correspond to a coordinate trans- 


formation. Instead it seems preferable to deal in terms of 


perturbations of the physically significant quantity, the curvature. 


2, Notation 
Space-time is represented as a four-dimensional Riemannian space 


with metric tensor of signature «2.  Covariant differentiation 


8 ab 
in this space is indicated by a semi-colon. Square brackets around 
indices indicate antisymmetrisation and round brackets symmetrisation. 


The conventions for the Riemann and Ricci tensors are:- 


Ve; (be) - a Brock Vs , 
Aat = ¿A 


D beca is the alternating tensor, 


| Units are such that k the gravitational constant and c, the speed of 


light are one. 


3, The Field Equations 


We assume the Einstein equations: 
Rab - = 3.4 R = ~“ lab 


where Tob is the energy momentum tensor of matter. We will assume 


that the matter consists of a perfect fluid. Then, 


Tab = fa U4U, 4 f hab 


where U, is the velecity of the fluid, U, i uj 
p is the density 
h is the pressure. 
hob = gut UL Uy is the projection operator 


into the hyperplane orthogonal to [UN 
NE so, | 
We decompose the gradient of the velocity vector [UN as 


where Ue = Usb u* is the acceleration, 

Bs °° is the expansion, 

Tar < (cie) ha ie - 5 hab is the shear, 

Wa = utd hz hi is the rotation of the 
flow lines Us: We define the rotation vector Wa as 

Wa * E Nabed itu 
We may decompose the Riemann tensor B bod into the Ricci tensor 
Rab and the Weyl tensor Cp! 


R abed = Cabcd ia Sata Re] T S bfe Raja = Ry Seats YA) b 
Cabed = Clab]Eda) 


) 
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C bcd is that part of the curvature that is not determined locally by 


the matter. It may thus be taken as representing the free gravit- 
ational field (Jordan, Ehlers and Kunat い 5)). We may decompose it 


into its "electric" and “ma netic" components. 
i Bg 


E.e = © Cobpq UPUT , 


= fs 5 
Hab üi * Curt Y qrbs Upu ? 
cd = tera [c d] 
C. < BUp Eq - 46 Ey 


- 2 Nabed u? Hale ui. E ve Hs ta uy 


Par = Eco » Wa Hey > 


E b and Hob each have five independent components. 


We regard the Bianchi identities, 


Li tedyei xdi. 


as field equations for the free gravitational field. 


Then 


C dia Pt = -R c[bia]+ t SepoR a] (Kundt and retimper, (64) 


analogous to the Maxwell equations. 
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Using the decompositions given above, we may write the 


ha” H be: d p? - 3 ab > m YJab col Qa ge. È に (K+ h) Wa 
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where | indicates projection by hy orthogonal to 
Gef. uta, MA, 


The contracted Bianchi identities give, 


(Rus Ry" a e Tola 
Ji + (ga )8 = © ; 
(+ f) ka + huso Pia = 0 ・ 


The definitien of the Riemann tensor is, 


U a [be] = 2 Rapbe uP = 


Using the decompositions as above we may obtain what may be r 


as "equations of motion", 


Ue 
a 


se in a form 


(1) 


(4) 


egarded 


B = 20 -20°-50 + uam = Peer , 


u 


l Sap E ab T Wac 


: 2 
aV OU. 053 FI 9 


\ a x . a 
a | Fii (a «x - 2g"4 tà 456) + い 。 い し 


. P 
+ Utp) ha 


n 2 a 
where 260"= Wab W 


We also obtain what may be regarded as equations of constraint. 


8 hi = 2 (( wre,b+ oou) hz - At (Wee + Cal] 
Da” = 2 GO a Ta > 

hi utc + の de | 
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We consider perturbations of 


ht 


2 b 
h ; 20 = Oab gs 
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in conformally flat, that is 


Cabcd 


By equations (1) - (3), this implies, 


も = Wab= 0 


b 
har» = 03 6.4, Na 


universe that in the undisturbed 


(11) 


(12) 


state 


If we assume an equation of state of the form, れ = Te G2 j 


gnen by (6), (10), typ hPa = Q0 = i è 


This implies that the universe is spatially homogeneous and isotropic 
since there is no direction defined in the 3-space orthogonal to U 

In this universe we consider small perturbations of the motion 
of the fluid and of the Weyl tensor. We neglect products of small 
quantities and perform derivatives with respect to the undisturbed 
metric. Since all the quantities we are interested in with the 
exception of the scalars, ps To 6 have unperturbed value zero, we 
avoid perturbations that merely represent coordinate transformation 
and have no physical significance. 


To the first order the equations (1) - (4) and (7) - (9) are 


Eu th Pjb > (13) 
Hoy? = (urh)wa , (41) 
En PEE + hc ede U pn -y (poh) Tab; (15) 
Bun HO bc ES wo , (16) 
O =- 78+ uo Fs) > (17) 
Coan 2 7$ + ALP Made, (18) 
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+ V PI 
Ga E Ey wh s ha + Uppa Nb . (19) 


from these we see. that perturbations of rotation or of E, or H do 


not produce perturbations of the expansion or the density, Nor do 


produce rotational perturhations. 


perturbations of Eb and Hob I 


| 4, The Undisturbed Metric 
Since in the unperturbed state the rotation and acceleration 

are Zero, Ua must be hypersurface orthogonal. 
Ua. * Tae 


a 


H 


where ‘© measures the proper time along the world lines. As the 
surfaces て = constant are homogeneous and isotropic they must be 
3-purfaces of constant curvature. Therefore the metric can be 


written, A 
ds* =-dr*+ Q dy” 


where Di Q (©) ; 


dy* is the line element of a space of 


zero or unit positive or negative curvature. 


We define t by, 


dr i 
2 
then dsta Q(-dt + dy?) 
In this metric, Uy = ( - 6, O, O) o) 


& @e Mbe 22 


y. 


(prime denotes differentiation with respect to t) 


Then, by (5); (7) 


b=- leh g a (20) 


sa = -7(4+3h) , (21) 


If we know the relation between u and PL y we may determine El 
We will consider the two extreme cases, fa = O (dust) and h= E 
(radiation). Any physical situation should lie between these, 


For fh = 0 


By (20), ^ “s En PA = const. 
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E represents the energy (kinetic + potential) per unit mass. 
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If it is non-negative the universe will expand indefinitely, other- 


Wise it will eventually contract again. 


By the Gauss Codazzi equati 


hypersurface て = cons 


d vc 
ak: 

に = の X 
T um. K 


Mm 
AÑ 
o 

È 

rs 


(2 = E sinh t 5 


b. 


By (6) 


Rotational Perturbati 


m|- 


n|- 


ons *R,the curvature of the 


Rea (ag の +) 


NA = =$ 
< a 
AD Am 8 
"I TA el e 


C 
> "A = 6 » 
/ * > © 
cos È - 1) << 
3 : es 
m _ Wot n 


Thus rotation dies away as the universe expands, This is in fact a 


Statement of the conservation of angular momentum in an expanding 


universe, 


6. Perturbations of Density 


For h = 19) we have the equations, 


p = - コン 


These involve no spatial derivatives. Thus the behaviour of one 
region is unaffected by the behaviour of another. Perturbations 
Will consist in some regions having slightly higher or lower values 
of E than the average. If the universe as wh^le has a value of E 


8reater than zero, a small perturbation will till have E greater 
? に 


than zero and will continue to expand, It will not contract to 


form a galaxy. If the universe has 


small perturbation can contract, 


contracting at a time 


contracting, where 
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T is the time at which the 
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There is only any real instability when 


measure zero relative to all the po 


However this cannot really be 
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as there might be some reason 


8 © earlier t 
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why the universe should have E 


of E less 
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a value 


zero, a 


However it will only begin 


han the whole universe begins 


universe begins contracting. 
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ible values E can have. 
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Thus the porturbation grows only as て 7 う , This is not fast enough 
LO produce galaxies from statistical fluctuations even if these 


could occur, However, 
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in the early stages This makes it 


Statistical fluctuations to occur 
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ps before, a perturbation cannot contract unless it has a negative 


value of E. The action of the pressure forces make it still more 


difficult for it to contract. Eliminating 0, 
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to our approximation, 


ec b こ ー 
h Vy haVb is the Laplacian in the hypersurface て = constant. 
E yo 
st ) of 
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fe represent the perturbation as a sum of eigenfunetions 


this operator, where, go © un 
sa th = © 
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These eigenfunctions will be hyperspherical and pseudohyperspherical 


harmonics in cases (c) and (a) respectively anc plane waves in case 


(b). In case (c) n will take only discrete ` lues but in (a) and 


(b) it will take all positive values. 


where Ho is the undisturbed density. 
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These perturbations grow for as long as light has not had time to 
travel a significant distance compared to the scale of the perturbation | 
( vw = i, Until that time pressure forces cannot act to even out 


perturbations. 
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Ye obtain sound waves whose amplitude decreases with time. These 
Zz 
results confirm those obtained by Lifshitz and haat ko , 
From the forgoing we see that galaxies cannot form as the result 


of the growth of small perturbations. We may >xpect that other non- 


gravitational force: will have an effect smaller than pressure equal 


to one third of the density and so will not cause relative perturbations 


to grow faster than T , To account for galaxies in an evolutionary 
universe we must assume tnere were finite, non-statistical, initial 


inhomogeneitics. 


7. The Steady-state Universe 


To obtain the steady-state universe we must add extra terms to 


(10) 


the energy-momentum tensor. Hoyle and Narlikar use, 
i . d 
| ga = HM UL Ub + h hu x Cu Cr * ノ 2 Gob Cdl ` (20) 


where, Ce & & 


T. - (p+ fr) U a 


i sb 
Since Tap’ -o 


ん Y (H h) 8 + ue Lia p. te =. 0 


(phi) eat Push ao ha Ci Cao, 


There is a difficulty here, if we require that the "C" field 


ghould not produce acceleration or, in other words, that the matter 


created should have the same velocity as the matter already in 


existence,We must thén have 


A Cr = 0 


However since C is a scalar, this implies that the rotation of the 

medium is zero. On the other hand if (23) does not hold, the equations 

are indeterminate (c.f. Raychaudhuri and Banner jee ). In order to 

have a determinate set of equations we will adopt (23) but drop the 

requirement that Ca is the gradient of a scalar. The condition (25) 
€ 

is not very satisfactory but it is difficult to think of one more 


satisfactory.  Hoyle and Narlikar (12) seek to avoid this difficulty 


by taking a perticle rather than a fluid picture. However this has a 


(13)y, 


serious drawback since it leads to infinite fields (Hawking 


From (17), k 


Ca ## Ba | da E+E s h) 6 


Thus, small perturbations of density die away. Moreover equation (18) 


still holds, and therefore rotational perturbations also die away. 


Equation (19) now becomes 


à T TENES ; ‘ 14 | 
These results confirm those obtained by Hoyle and Naka B. We 
see therefore that gslaxies cannot be formed in the steady-state 
universe by the growth of small perturbations, However this does not 


d- Y 


exclude the possibility that there might by a self-perpetuating 


system of finite perturbations which could produce galaxies. 
(15) 


e A, 7 TD rh nr コ c ALMA (1 6) 
(Sciama Roxburgh and Saffman ) 


8. 


4 


We now consider perturbations of the Weyl tensor that do not 


1 


arise from rotational or density perturbations, that is, 


tb b 
E ; i 
E ab = | | a h o 
c &77 P E et ybs M 
Multiplying (15) by U Ve, and (16) by hr la Ue Me 5 


we obtain, after a lot of reduction, 


Eno - (Eede NY n hti): HH + E ab O 
+ Eab( Ò t £0 + (psh)) + Tas [ 5 Gu fo) MENTE fa] -—- fon’ 


A è a è 
In empty space with a non-expanding congruence U^ this reduces to 


the usual form of the linearised theory, 


2, 


[I Es <6 


The second term in (2h) is the Laplacian in the hypersurface 


^t = constant, acting on Eb 。 We will write E ab as a sum of 
eigensfunctions of this operator. 
- - n) (m) 
E b 2 ZL A V. b 
o. 
e i, 
where Vab =0 zy 
j ed e bis NM (n) 
A la ¢ hs hi) Whe? h b. s pud 2. ab m" 
? (2 
V sb a 
ab =D MV P = 0 
ES こい In) (m) 
Then Esp = 2 らき Mau. 
n HD + 
= i Kn) m) "Y í an 
E ab の 〆 | SW AE ^ M 
L Cè w LAN 


N) n) 
gimilarly, Ca C D Mi 


when by (19) = 


Substituting in (21) 


; ve " T " " 12 
n" y os AU, ae? |^ + 30 e 6 de. (Xx 3p) 2] 


+ p? (Quad) ++ Qt f) = © 


We may differentiate again and substitute for D', 
For n! 


= zi... 
and Cl. q 


T 
AU ge Po. yA 


80 the gravitational field Bap decreases as Q hs 


and the "energy" 


ab ab 
BU +8, H ) 


Bianchi identities may be written, to the linear approximation, 


mn -6 : -— : 
as Q 。 We might expect this as the 


n PA (qc abit.) dala 


Therefore if the interaction with the matter could be neglected 


C A D xt | } anó 7 ー1 
E od would be proportional to C) and E b ; Hop to (1 o 
In the steady-state universe when pH and 6 have reached their 
113 j alea TA = E r f 
equilibrium values, Ra. = [5 t ^j 


Jae Recta; € Seta Rsb] 


z: 0 


mhus the interaction of the "C" field with gravitational radiation is 


equal and opposite to that of the matter. There is then no net 


; " ‘ ー1 
interaction, and Bap and Hb decrease as GL 


1 a 
ife we o dab 
2 (E È "ab i Tab ) 


of the metrie. It is therefore proportional to the frequeney squared 


The "energy" depends on seeond derivatives 
times the energy as measured by the energy momentum pseudo-tensor, in 
a local co-moving Cartesian coordinate system which depends only on 
first derivatives. Since the frequency will he inversely proportional 
to C2 ; the energy measured by ihe pseudo-tensor will be proportional 


to C 4 as for other rest mass zero fields, 


9, Absorption of Gravitational Waves 
As we have seen, gravitational waves are not absorbed by a 
perfect fluid. Suppose however there is a small amount of viscosity. 
We may represent this by the addition of a term A Tab to the 
energy-momentum tensor, where "A is the coefficient of viscosity 

(17)), 


(Ehlers, 


Since Tal SB 


We have 


ptr )e-2»v€* zo (25) 


" ) m 
(uh) UU + ha, Aa PAGG" ho. = 


(26) 


gquations (15) (16) become 


し 


di 
M + Eabb + ha ede ut Hi, i = - E (i fo) Sab 


iù A ( Ebay 9) A (27) 
n {: x d:e 
Am | HO - h (abede us - =-+ Hap | (28) 


The extra terms on the right of equations (27), (28) are similar to 

conduction terms in Maxwell's equations and will cause the wave to 
T . . 

decrease by a factor e を と 。 Neglecting expansion for the moment, 


suppose we have a wave of the form, 
= x と や で 
に ab” E b € 


: : ; = ‘ 2- R 
This will be absorbed in a characteristic time A independent of 
frequency. By (25) the rate of gain of rest mass energy of the 


— 2 - 2 
matter will be 2540 ^ which by (19) will be 2X E V . Thus the 


-2 -2 

available energy in the wave is 4 oE VO, This confirms that the 
density of available energy of gravitational radiation will decrease 

0) -l |; T ; 1 
8S in an expanding universe. From this we see that 
gravitational radiation behaves in much the same way as other 
radiation fields. In the early stages of an evolutionary universe 
When the temperature was very high we might expect an equilibrium to 


bo set up between black-body electromagnetic ; diation and black-body 


gravitational radition. Since they both have two polarisations their 


energy densities should be equal. As the universe expanded they woul. 


both cool adiabatically at the same rate. AS we know the temperature 
of black-body extragalactic eleciromagnetic radiation is less than 
59x , the temperature of the blacx-body gravitational radiation musi 
be also less than this which would be absolutely undetectable. Now 
the energy cf gravitational radistion does not contribute to the 
ordinary energy momentum tensor Tab 。 Nevertheless it will have an 
active gravitational effect. By the expansion equation, 


4 KS 2. 
= ! 1 ^ 
0 =- z0 -20 -A(A95b) 
For incoherent gravitational radiation at frequency v , 
2 
i bid 
res: QE v 
But the energy density of the radiation is 4 Ev 


n 007-5$0-£Hs-(ursh) 


whei? Lu 


u, is the gravitational "energy" density. Thus gravitational 
G < | 5 n 


Padiation has an active attractive gravitational effect, It is 
interesting that this seems to be just half that of electromagnetic 
radiation. 


(18) 


3 F) 
and Hoyle and Harlikap 19), 


It has been suggested by Hogarth 
that there may be a connection between the absorption of radiation 
and the Arrow of Time. Thus in universes like the steady-state, in 


Which all electromagnetic radiation emitted is :ventually absorbed by 


Other matter, the Absorber theory would predic retarded solutions of 


the Maxwell equations while in evolutionary universes in which 


electromagnetic radiation is not completely absorbed it would predict 
advanced solutions. Similarly, if one accepted this theory, one would 
expect retarded solutions of the Einstein equations if and only if all 
gravitational radiation emitted is eventually absorbed by other matter. 
Clearly this is so for the steady-state universe since ^ will be 
constant. In evolutionary universes ^ will be a function of time. 

We will obtain complete absorption if [Ade diverges. Now fr? a gas, 
hoc r where T is the temperature. For a monatomic gas, peo, 
therefore the integral will diverge (just). However the expression 

used for viscosity assumed that the mean free path of the atoms was 
small compared to the scale of the disturbance. Since the mean free 
path x ule £2 73 and the wavelength & (7 , the mean free path will 
eventually be greater than the wavelength and so the effective viscosity 
will decrease more rapidly than (? -1 . Thus there will not be complote 
absorption and the theory would not predict retarded solutions. 


However this is slightly academic since gravitational radiation has not 


yet been detected, let alone investigated to see whether it correspond s 


to a retarded or advanced solution. 
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CHAPTER 5 


Gravitational Qadiation In An 


Expanding Universe 


Pad 


Gravitational radiation in empty asymptotically flat 
space has been examined by means of asymptotic expansions 
AAN " -—————— 
by a number of authors. lhey find that the different 
components of the outgoing radiation field "peel off", thet 
is, they go as different powers of the affine radial distance. 
If one wishes to investigate bow this behaviour is modified 


5 PIS 
E 


by the presence of matter, one is faced with a difficulty 
that does not arise in the case of, say, electromagnetic 
radiation in matter. For this one can consider the radiation 
travelling through an infinite uniform medium that is static 
apart from the disturbance created by the radiation. In the 


case of gravitational radistion this is not possible. tor, 


would cause it to contract in on itself and it would cease to 
be static. Hence one is forced to investigate gravitational 
raGiation in matter that is either contracting or expanding. 
äs in chapter 2, we identify the Weyl or conformal 
tensor ES: with the free gravitational field and the 


() 


ae . » . m .. r è a 
Ricci-tensor hab with the contribution of the matter to the 


curvature. Instead of considering gravitational radiation in 


asymptotically flat space, tnat is, space that approaches 


flat space at large radial distances, we consider it in 
asymptotically conformally flat space. As it is only 
conformally flat, the Ricci-tensor and the density of metter 
need not be zero. 


To avoià essentially non-gravitational phenomena such 


as sound waves, we will consider gravitationalradiation 


Aw 


travelling through dust. It was shown in Chapter 2 that 

conformally flat universe filled with dust must have one of 
the metrics: 
> 2 Oca Bu etr 

(n) det, QM db du rn p LU |n OF y) 

Sh Bee c) TES 
| ae a 2 THEM LIA 
A nar E Lp = a (LE + fia “Od q )| 
( へ 2 pa / 

bs are Ch 2) 


à fa? c? 


cl. 2 $) : dt i B dud eink / CLO rS V Kil // 
SL < A (cesh.& - 1) (1.9) 


Type(a) represents a universe in which the matter 


expands from the initial singularity with insufficient e 


to reach infinity and so falls back again to another 


singularity. It is therefore unsuitable for a discussion of 


gravitational radiation by a method of asymototic expansions 


since one cannot get an infinite distance from thie source. 

Type (b) is the Zinstein-2e Sitter universe in which the 
matter nas just sufficient energy. to reach infinity. It is 
thus a special case. D. Norman (2) has investigated the 
"peeling off" behaviour in this case using Penrose's conformal 
technique (6). He was however forced to make certain assumpt- 
ious about the movement of the matter which will be shown to 
be false. Moreover, he was misled by the special nature o: 
the Hinstein-De Sitter universe in which affine and luminosity 
distances differ. Another reason for not considering radiation 
in the tinstein-Ve Sitter universe is thet it is unstable. 


the passage of a gravitational wave will cause it to contract 


asain eventually and develop a singularity. 


4e will therefore consider radiation in 8 universe of 
type (c) which corresponds to the general case where the 


matter is expanding with more than enough energy to avoid 


contracting again. 


2. The Newman-Penrose Formalism 
(3) 4 


we employ the notation of Newman and Penrose. 


た mi^ — 
tetrad of null vectors,U, MU, M m^ 


is introduced 


gravitational radiation by a method of asymototic expansions 


since one cannot Set an infinite distance from thi@ source. 
Type (b) is the Uinstein-De Sitter universe in which the 
matter nas just sufficient energy. to reach infinity. It is 


EE 
thus a special case. D. Norman (2) has investigated the 
Penrose's conformal 


(6 $ a ; 
‘ è, He was however forced to make certain assumpt- 


technique 
ions about the movement of the matter which will be shown to 
be false. Moreover, he was misled by the special nature of 
the Kinstein-De Sitter universe in which affine and luminosity 
distances differ. Another reason for not considering radiation 
in the winstein-De Sitter universe is that it is unstable. 
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the passage of a gravitational wave will cause it to contract 


azain eventually and develop a sinzularity. 

ie will therefore consider radiation in a universe of 
type (c) which corresponds to the general case where the 
matter is expanding with more then enough energy to avoid 


contracting again. 


2. The Newman-Penrose Formalism 


we employ the notation of Newman and Penrose, 63) A 


[^ mi^ E 
tetrad of null vectors ,U, N, M, m^ is introduced 


AM. y I g 
where: - A, a t PT mi- し 。 fm! a - È 


we label these vectors with e tetrad index 


z^ . (Ln m^ M" ) 421,2,3,4 
tetrad indices are raised ana lowered with the metric 
ab 23 - 
/ : ) : o i © GO 
ob | 
[fo o o (2.1) 
oo 0 «| 
O Q | O 
L ] 
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we have C PT " y 2” 3" 
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MV pv ù ¿Y P ur £z 
s b Ek D m ie PM (23) 


Ricci rotation coefficients are definod by: 


bc "m 

Y = Zpw 2! 2° (i.d) 
a 
abe boc 


complex combinations of rotation coefficients defined 8S 
follows: | 
K - Y = Ly; MEL” 
31 
= ^P 
Tt = 23 z fiy. T m E 
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C L ( A Y) s Lage b 8] 
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In fact it is more convenient to work in terms of twelve 


p- [A 


(s (Eo) e MAR poli 


/ 


/ 
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V 


\ 
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4. Coordinates 


Like Newman and Penrose, we introduce a null coordinate 


(= K") 
( I^ V ; a 
q Uy wiv の tE) 


vV 
we take / = Ul e Ad 
press pe / / 


geodesic and irrotational. This implies 


Ko 


hus た ん will be 


qa - E = A t p - 3 と ex 7 11-4 ーー D T 
we take n^, mé [^ ^ to be parallelly transported 
/ 


along た 。 This gives , 
m ec A 
Tt 3 し zx C la 
AS a second coordinate we take an affine parameter ràe x” 


long the geodesics L 


v, LP 3 | G y) 
3 P4 


X and XY aro two coordinates that label the geodesic 


in the surface U = const. 


O (3.5) 
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In these coordinates > bin < $ AA b di. 


j^ a E ux i 1 
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pxpressing the rotation coefficients in terms of the metric, 


we have: 


) a ; E E ¢ と (s 3) 
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As in Chapter 2 we use the Bianchi identities as 


field equations for the Weyl tensor. In the Newman-Penrose 
formalism they may be written: 
(I am indebted to R. G. McLena;han for these) 
9 Y, -Dy + Don ~ 29 = qa Ye ー abd Cla qe) の 。 
map Pot 25 de 6.31) 
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4, The Undisturbed Metric 


The undisturbed metric may be written 
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(ls? = sL ae” - q, f 2- sinhto (LOs in*d p> | 
の = A€cosh € -1) 
put la = ら -ø 
d 9? s ^ | -du T Idudt - Sin | Gies Nd Ervin ed // 
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chen 


ù is a null coordinate 
To calculate f“, the affine parameter, we note thit C 

is an affine parameter for the metric within the square 
brackets. Therefore T c j “al, Cu 6 Cu, e の ) (4 2) 
will be an affine parameter for (L*1) 

B is constant along the null geodesic. Normally it 
would be taken so that | - O when に ad. 。 However, 
in our case it will be more convenient to make it zero and 


+ = RE (4 3) 


define l^ as 
This means that surfaces of constant Y are surfaces of 
constant と . This may seem rather odd, but it should 

be pointed out that the choice of R will not affect the 


symptotic dependence of quantities. That is, if 


" {= O(v7^) 
(reped) pis rob 


J 


Then 


It proves easier to perform tre calculations with this 


choice of た but all results could be transformed back 


to a more normal coordinate system. 


Prom (4.3) T = R^ fE sinh) =I sinh tr i£] G 4) 


The matter in the universe is assumed to be dust so its 
energy tensor may be written 
ご Va Y (4. S) 
T ub P C b 


For the undisturbed case, from Chapter 2 
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Therefore if we try to expand p as a series in power of © 
the result will be very messy and will involve terms of 


Y) 
the form loq S * 


*[t should be pointed out that the expansions used will 
only be assumed to be valid asymptotically. They will not 
be assumed to converge at finite distances nor will the 
quantities concerned be assumed analytic. (see A. Erdelyi: 


Asymptotic Expansions - Dover 


This does not invalidate it as an asymptotic expansion but 


it makes it tedious to handle. For convenience therefore, 


we will perform the expansions in terms of Jot) which 
will be defined in general as the same function of の as 
it is in the undisturbed case. ‘hat is | 


S « D (cosà E-i) ) 
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then pe 
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For the third and fourth coordinates it is more convenient to 
use stereographic coordinates than spherical polars. 

Since the matter is dust its energy-momentum tensor 
and hence the Ricci-tensor have only four independent 
components. We will take these as Ad. Be e 
(since Po, is complex it represents two components) 
In terms of these the other components of the Ricci-tensor 


may be expressed as: 
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For the undisturbed universe with the coordinate system 


given: A - E " A . 
LY (52.7 
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Using these values and the fact that in the undisturbed 


è S 
universe all the Kd are zero, we may integrate equations 


(3. 10-50) to find the values of the spin coefficients for the 


unperturbed universe: Q Q (4 L et 
si i 


ELE € AA WE eun ' 
[2 -e**)5x5 A(E- Jo iat d 

Ew Ga t= v S x X Du 

gie del di, a(S Be INS. 


| 


La い / peo 
ER =P S Se 一 7AVSe" 


T 


f^ 


EE 
> 
t 
| 
JE 
し 
t 
(で 
N 
と 
a 
Es 
いい 
E 
js 
C^ 
= 
+ 
Ko 
G 
C 
J 


| 
U wo d we a 


— 


5. Boundary Conditions 


Ye wish to consider radiation in a universe that 


asymptotically approaches the undisturbed universe given 
above. (b and /\ will then have the values 
OD 

given above plus terms of smaller order. To determine this 
order and the order of (D. and % , tuere-are two 
ways in which we may proceed. ie may take the smallest orders 
Ghat will permit radiation, that is A = Ol) e 
Larger order terms than these in D, 9 /\ and 

D turn out to have their VA derivatives 
IMS only on themselves and not on the r^ coefficient 
of Y , the radiation field. They are thus disturbances 


not produced by the radiation field and will not be considered. 


Alternatively we may proceed by a method of successive 


£ 


approximations. We take the undisturbed values of the spin 


coefficients and use them to solve the Bianchi Identities ^s 
field equations for the conformal tensor using the flat space 
boundary coridition that Ve + olr”) . Then substituting 
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hold on one hypersurface (ù = const.) they will hold on 


succeeding hypersurfaces and that these conditions are the 


most severe to permit radiation. 
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and not on the radiation field. it tnerefore represents a 
type of disturbance unconnected with radiation. If it is 


zero on one hypersurface, it will remain zero. In this case 
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The asymototic group is the group of coordinate transformations 
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Sinzularities 


If the binstein equations without cosmological constant 
are satisfied, a Robertson-Walker model can 'bounce' or avoid 
a singularity only if the pressure is less than minus one- 
third the density. This is clearly not a property oossessed 
by normal matter though it might be possessed by a field of 
negative energy density like the 'C' field. However there is 
a prave quantum-mechanical difficulty associated with the 
existence of negative energy density, for there would be 
nothin, to prevent the creation, in a given volume of space- 
time, of an infinite number of quanta of the negative energy 
field and a corresponding infinity of particles of positive 
energy. If we therefore exclude such fields, all Robertson- 
Walker models must be of the 'big-bang' type, that is they 
have a singularity in the past and maybe one in the future 
as well. It has been suggested 1 that tne occurrence of 
these singularities is a consequence of the high degree of 


symmetry of the Robertson-Walker models which restricts the 


expansion and contraction so that they are purely radial and 


that more realistic models with fewer or no exact symmetries 


would not have a singularity. This chapter will be devoted 


to an examination of this question and it will be shown that 


provided certain physically reasonable conditions hold, any 
model must have a singularity, that is, it cannot be a 
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This also shows that singular points of congruences are 


points where neighbouring geodesics intersect. 
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is the energy density in the rest-frame of the observer and 
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Condition (a) will be satisfied by a perfect fluid with 
oper 

density M 2O and pressure P p zp 。 It implies kV V70 for 
any time-like or null vector i A Therefore by equations 
(1) and (5) any time-like or null irrotational geodesic 
congruence must have a singular point on each geodesic withi 
a finite affine distance. Obviously if the flow-lines form 
on irrotational geodesic congruence, there will be a physical 
singularity at the singular points of the congruence where 
the density and hence the curvature are infinite. ‘this will 


be the case if the universe is filled with non-rotating dust 


However, if the flow-lines are not geodesic (ie. non-vanishing 
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pressure gradient) or are rotating, equation (1) cannot be 
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The Robertson-welker models are spatially homogeneous 
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and isotropic, that is, they have a Six parameter group of 


motions transitive on a spece-like surface. If we reduce the 
symmetry by considering models thet are spatially ^omozereous 


but anisotropic (thatiis, they have a three parameter group 


of motions transitive on a space-like hypersurface) then 


the matter flow may have rotation, acceleration and shear. 


Thus there would seem to be the possibility of non-singular 
4 ee ee eee + i La を 

models. L. Shepley has investigated one particular | 

homogeneous model containing rotating dust and has shown that 


there is always a singularity. Here a general result will be 


proved. 
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There must be 8 Singularity in every 
condition (a) and, 
(b) there exists a Gy of motions on the space or on 


universal covering space , -- which is transitive on at 


t vee section 5 


1 bats の マテ 


least one space-like surface but space-time is not stationary, 


(c) the energy-momentum tensor is that of & perfect fluid, 
a 
y ot = (jv P) u U > PO ob è U is the vanzent to 


flow-lines and is uniquely defined as the time-like eigen- 


vector of the kicci tensor. 


R, the curvature scalar must be constant 


D 
Ta ar 


like surface of transitivity リア of 


B * 4 e . $ . . * っ G 
nust be in the direction of the unit time-like normal V 


vo. A” . 


wnere E Re z f »o 


CR. is an indicator 
C CS sa. 


Il 


+1 if JT is past directed 


~ 


Il 


-1 if Ria is future directed 


Then Va ; b1 =0 e 
Thus Va is à congruence of seodesic irrotational time-lilce 


a b 
vectors. By condition (a), ea Vv >O à 


therefore the congruence must have a singular point on each 
slo ee de 


geodesic ( by equation 1) either in the future or in tne vast. 


jurther, by the homogeneity, the distance along eacn 


3 
from dl” to the singular point must be the same. for each <eo.esic. 


thus if the surfaces of transitivity remain space-like, 


: 4 } n "o, NE. 1 4 
must degenerate into, at the most, a 2-surface C^ which will 


be uniquely defined . Let M be the subset of the flow-lines 


2. 
of the matter which intersect E » Let し be the non- 


。 wince there is a 


: 3 
empty subset of H intersected by M 


p ^ 


3 
flou-lines trough li must intersect the 2-surface C^. ‘thus 


the density will be infinite there an: there will be a 
physical singularity. Alternatively if the surfaces of 
transitivity do not remain space-like, there must be at 
least one surface which is null - call this S^. At ^ FE = 0。 
Rio. 1 O ( irk. is Zero, we can take any other scalor 
oolynomial in the curvature tensor and its covariant derivatives. 
iney cannot all be zero if space-time is not stationery). we 


introduce a geodesic irrotational null congruence on S^ with 


Y / 


tangent vector go where b. a Then by eouation (5), 
there will be a singular point of each nuli geodesic in 
within finite affine distance either inthe future or in the 
past. “he 2-surface of these singular points will be uniquely 
defined. The same argument used before shows that the density 
becomes infinite and there is a physical singularity. in fact 


as S? is a surface of homogeneity, the whole of S^ will be 
Singular and it is not meaningful'ito cali it null or to 
Gistinguish tnis case from the case where the surfaces of 
transitivity remain svace-like. 


The conditions (a), (b), (c) may be weakened in two ways. 


Condition (b) that there is a group of motions throughout 


; à / i XA 
space-time may be replaced by (b') aná (d). 


/ 


TE 
(D) Thee is a space-like hypersurface H in which there 


| 3539 T EE PE a i ^ | ES dl Ps a A yo a } dis zr" 
are three independent vector fields X such that 


e ^ Q 
= O . E 3 
d des Y^ bede ^ O on H 。 That is, there 
i ^ m x | 
| i$ one homogeneous space section. 
(à) ‘There exist equations of state such that tho vauchy | 
À = ME 。 " "TEN 
development of is determinate. 
Then succeeding space-like surfaces of constant R 
are homogeneous and much the same proof can be given that there 
/ 
ix vis nino PON | des a pA > 5 (rh) fan) d. 
are no non-singular models satisfying (a), (b), (oc), (a). 
The only property of perfect fluids that has been used 
in the above proof is that they have well defined flow-linos 
intersection of which implies a physical singularity. Obviously) 


however, tais property will be possessed by a much more gen 


class of fluids. For these, we define the flow vector as 


time-like eigenvector (assumed unique) of the energy-momentum 
tensor. ‘men we can replace condition (c) on the nature of 


the matter by the much weaker condition (e). 

(e) If the model is singularity-free, the flow-lines form 
à smooth time-like congruence with no singular points with 
a line tarough each point of space-time. 


Condition (e) will be satisfed automatically if conditions 


(a) and (c) are. 


This proof rests strongly on the assumption of 


homogeneity which is clearly not satisfied by the physical 


universe locally though it may hold on a large enough scale. 
However it would seem to indicate that large scale effects 
like rotation cannot prevent the singularity. 

it is of interest to examine the nature of the singularit 
in the homogeneous anisotropic models since this is more 
likely to be representative of tue general case than that of 
the isotropic models. It seems that in general the collapse 
will be in one direction,” that is, the universe will collapse 
down to a 2-surface. Near the singularity, the volume will be 


proportional to the time from the singularity irrespective of 


the precise nature of the matter. It also appears that the 


テン 


nature of the particle horizon is different. There will be 
a particle horizon in every direction except that in which 


the collapse is taking place. 


4, Singularities in Inhomogeneous Models 


o Cm - : 6 
Lifshitz and Xhalatnikov 


claim to nave proved that & 
general solution of the field equations will not have a 
singularity. ‘Their method is to contract a solution with a 


singularity which they claim is representative of the 
o - 


general solution with a singularity, and then show that it 


has one fewer arbitrary function than a fully general solution. 


Clearly their whole proof rests on whether their solution 


they give no proof. 


is fully representative and of that 
indeed it would seem that it is not representative since it 
involves collapse in two directions to a 1-surface whereas 
in general one would expect collapse in one direction 
2-surface. In fact their claim nas been proved false by 
Penrose’ for the case of a collansing star using the notion 
of a 'closed trapped surface'. A similar method will be 
used to prove the occurrence of singularities in 'open ' 


universe models. 
5. 'Open' and 'Closed' Models 


The method used by Penrose to prove the occurrence 
of a physical singularity depends on the existence of a 
non-compact Cauchy surface. A Cauchy surface will be taken 
to mean a complete, connected egpece-=like surface that 
intersects every time-like and null line once and once only. 


Not all spaces possess a Cauchy surface: examples of those 


1 3 i s : の 1 n 2 " = 9 
that do not include the plane-wave metrics, the Godel model,” 
p~ ー ュ マ m co 10 y E p do e p を < ^ remi ne 
and N.U.T. space. However none of these have any physical 
Significance. Indeed it would seem reasonable to demand of 
any physically realistic model that it possess a Cauchy 


surface. If the Cauchy surface is compact, the model is | 


commonly said to be 'closed'; if non-compact, it is said to 


to be 'open'. The surfaces, t = constant, in the Robertson- 


Walker solutions for normal matter are examples of Cauchy 
surfaces. If k = -1, they have negative curvature and it is 
frequently stated that they are non-compact. ‘nis is not 

| necessarily so: there exist possible topologies for which 
they are compact. However, the following statements may be 
made about the topology of the surfaces t = constant. 


If the curvature is negative, K = -1, the universal 


13 


SS 


covering Space is non-compact and is diffeomorphic to E 
Any other topology can be obtained by identification of 
points. ‘thus any other topology will not be simply connected 
and, if compact, must have elements of infinite order in the 
fundamental group. Further if compact, they can have no 
group of notions. 

If the curvature is zero, K = 0, tne universal covering 
Space is à?. There are eighteen possible topologies. 2 If 
compact they have a Gz of motions and Betti numbers, B, = S. 
"UI E 

If the curvature is positive, K = +1, the universal 


D 


covering space is S?. Thus all topologies are compact. ‘the 
Be 


do 4^ ^ ^" "n 42 
tti numbers are all zero. 


Since a singularity in the universal covering space 


implies a singularity in the space covered, Penrose's method 


is applicable not only to spaces that have a non-compact Gauci 


surface but also to spaces whose universal covering space 


has a non-compact Cauchy surface. ‘hus it is applicable to 


models which, at the present time, are homogeneous and iso- 


tropic on a large scale with surfaces of approximate homo- 


geneity which have negative or zero curvature. 


" 
Let ?^ be a 3-ball of coordinate radius r in a 3-surface 


2 (t = const.) in a Robertson-Walker metric with K = 0 or -1. 


DI 


a i , - 2 3 
Let q” be the outward directed unit normal to T , the boundary 


ii 


Te ME ; a ^ bi T i " 
of 9, in y? and let V" be the past directed unit normal to 


u^. Consider the outgoing family of null geodesics which 


; a cx ci am T 
intersect T" orthogonally. at T, 。 their convergence will 


e por Au t uU XS ehe EMEP) 


; A ; a TES v 
where SA £t are unit space-like vectors in i^ orthogonal 


: a i á 
to q and to each other, 


l 
2 [dh =k ~ + iw | 


ET 


t3 


terefore f 


li 


LE ja? O and K O or -1, by taking r large enough, we may 
sand d negative at T^. Therefore, in the language of Penrose, | 


2 
¡€ 


PY” is a closed trapped surface. 


Another way of seeing this is to consider the diagram 


in which the flow-lines are drawn at their proper spatial 


distance from anobserver. They all meet in the singularity 
t O 


at t = 0. If the past light cone of the observer is drawn 


on tais diagram, it initially diverges from his world-line 


( f «OQ ). It reaches a maximum proper radius ( f sà > 
and then converges again to the singularity ( po ). The 


intersection of the converging light cone aná tne surface H 


< 42 -— " rain " 
gives a closed trapped surface T^. If the red-shift of the 


cuasi-stellar 2C9 is cosmological then it will be beyond the 
point d è O if we are living in a Robertson-valker type 
universe with normal matter. However, the assumptions of 
homogeneity and isotropy in the large seem to hold out to tne 
distance of 509. Thus there is sood reason to believe that 
our universe does in fact contain a closed trapped surface. 
it snould be pointed out that the possession of a closed 
trapped surface is a large scale property that does not depend 
on the exact local metric. Thus a model thet had local irregul- 
arities, rotation and shear but was Similar on a large scale 
at the present time to a Robertson-Walker model would nave a 
closed trapped surface. 

Following Penrose it will be shown that space-time has 


a singularity if there is a closed trapped surface and : 


(f) E Y O for any observer with velocity 


(g) there is a global time orientation 


& non-compact Cauchy 


(h) the universal covering space has & 
surface H^. 
PROOF 
Li 
4 pie A dé ais a ST a 4 
Assume Space-time is singularity free. Let ? be the 
— i A n の “a 
set of points to the past of H^ that can be joined by a 
2 . . t . 
or its interior 


C 09 
SHOW 


smooth future directed time-like line to T 
RS "UE E. è n4 STAEN a 22d 
T Let B^ be the boundary of F 。 Local consderations 
that B^? - $^ is null where it is non-singular and is 
generated by the outgoing family of past directed null geodesic 
segments which have future end-point on T" and past end-point 

where or before a singular point of the null geodesic 


congruence.  oince at me , し ne convergence, p? 〇 and 
ayb 
since R pÉ C ヶ の by (f), the convergence must 
[^ も 


become infinite within finite affine distance. ‘thus B^ - y 
will be compact being generated by a compact family of compact 


ズ 
Hence B^ will be compact. 
5 ; : 1 
B^ arbitrarily closely 


then as follows: approximate B 
5 


Penrose's method is 


segments. 
smooth space-like surface and project B^ onto H^ by the 
normals to this surface. “his gives a many-one continuous 
a a 5 * i m dictom, o Zx 
mapping of B- into H^. since B^ is compact, its image RB 


ズ 
must be compact. Let d.(Q) be the number of points of B^ 
ly at the 


H^. UO) will change only 


mapped to a point Q of H 
ul ; eds 2 
intersection of caustics of the normals witn H^. Moreover, 


cL (0) can only change by an even number. 


by continuity 


[ 


3% 
Bus dir- 4 since this is the identity ana (HS た 


ana 


This is a contradiction, thus the assumption that space-time 


D 


is non-singular must be false. An alternative proceedure 


which avoids the slightly questionable step of approximating 


Z 
で 
Ly 


Y 


by a space-like surface is possible if we adopt condition 


WN 


(e) on the nature of the matter, then B^ may be projected 


continuously one-to-one onto u^? by the flow-lines.  -his 


i às NE E ; "RN. 
again Leads to a contradiction since p? is compact and H 


yn っ 1 3 om n3 > " 55 
-moophismn-—of—ly---into-—H-. 


ズ 
dizection away from H^ since ia tho direction 5owards-H^ 


8. singularities in 'Closed' Universes 


There is a singularity in every model which satisfies 
(a), (g) and (i). 
(i) There exists a compact Cauchy surface H 
o 

normal V has positive expansion everywhere on H^. 
PROOF 

yor the proof it is necessary to establish a couple of 
lemmas. Assume that space-time is singularity-free. The 
following result is quoted without proof, it mey irescily be 
derived from lemmas proved in reference 11. 

LI P and q. are conjugate points along a geodesic y 


and X is a point on y not in Po thenX must have a conjugate 


point in PY 


An immediate corrollary is tnat if q is the first point 


along y conjugate to p and y is in pa then y has no 

conjugate points in pq è. also since tne result that x 

has a conjugate point in pq can only depend on the valucs 

of A in pq, any irrotational geodesic congruence including 
MA 

the geodesic Y must have a singular point on y in pa. 

Thus if q is a point on mó and Y is the geodesic normal to 

M? through q, then a point conjugate to q along d cannot occur 

until after a point conjugate to M^. 

If M? is a complete connected space-like surface which 

intersects every time-like and null line from a point p, we 

may define a function over M? as the square of ihe geodesic 

distance from p which is taken as positive if the geodesic is 

time-like und negative if the geodesic is space-like. 'e call 

this the world function O with respect to p. For the closed 

set of values ぐ ZO » の will be a continuous ( in 

general multi-valued) function over w^. A time-like geouesic 
Ü from p will be said to be critical if it corresponds 

to à volue of Ó for which o ela O (i s i, 2, 3) 


4 
2 し 


一 


M i 
where (© are three independent vectors in m3, Clesriy 
し 


à critical geodesic must be orthogonal to M^. A geodesic 


which is critical will be said to be maximal if it corresponds to a 


local maximum of . 
Lemma 1. 
panta K ì vida or 2d 
A geodesic D cannot be maximal for a smooth M^ if there is a 


point X conjugate to w^ but no point conjugate to q on Y in qp, 
where q is the intersection of ^ and M^, 


. . / 
which vanish at X 


Let f and g be the Jacobi fields along 


m m 
and respectively. They may be written 
| 
n 
f = As) f/a, 
m mn 
n 
g = B(s)ga. 
m mn 
Then (AN : uml / 。 f must be positive for 
n " ( v a ノ n 
any h since if it were negative for any h by taking a =-bhh 


mb mn 
beyond q, it would be possible to have a point y on の beyond q 


X 


conjugate to before a point conjugate to P . If it were zero 


would be conjugate to P 。 This shows that the surface at q 


of constant geodesic distance from (PP lies nearer to fpf in every 


direction than the surface o£ の of constant geodesic distance from 


T à ^ Sa 
X does. Since > is conjugate to M^ the surface at q of 
f: 


constant geodesic distance from Vv lies closer to # in some 


direction than MÍ does, Hence 0 is not maximal. 


If M^ is compact or if the intersection of all time-like 
and null lines with M^ is compact, の must have a maximu.. 
A 
value, thus there must be a geodesic normal to M? through p 


longer than T 。 We use this to prove another lemis. 
Lemma 2 
If p lies to tne future (past) on a time-like geodesic 
U through gq, beyond a point z conjugate to a,and tnere 
" 


exists a compact Cauchy surface i^ through q, then there 

must be another time-like geodesic from p to q longer than y 
Let y be the last point conjugate to q on T before p. 

Let x be the nearest point to p conjugate to p in pq. Let 

r be à point in yx. Let KÍ be the set of points which have 


a future (past) directed geodesic of length rq from q. Then 


PR : . és | : n : 4. 
will be a space-like hypersurface through r. Let i 


be the set of points which nave at least one future (past) 


directed geodesic from q of length greater than rq. ‘hen the 


. i a ant T M AU ‘ 
boundary of F , Po K^. Since p is in F and since every 
past (future) directed time-like and null line from p intersects 


"C ; i " 4 à 4 T ds 
H^ which is not in F » they must also intersect 7^. Let L 


5 


be the intersection of J^ and these lines. Since iv” is 


compact L^ must be compact. Consider the function O° with 
È 3 p 


5 


respect to p over X^. Its maximum must lie in the compact 
region $e, But, by the vrevious lemma Y is not maximal, 
woreover, local considerations show that a singular point in 
the surface JÍ cannot be a maximum of の . ‘hus the 
1 maximum value of Q^ must occur for a geodesic from p ortao- 

gonal to b*. This must also be a goedesic from o to a of 
length greater than Y A 

Using these two lemmas tae theorem may be oroved. since 


H? are converging 


15 


everywhere on HÍ, there must be a point conjugate to H^ a 


the future (past) directed normals to 


finite distance along each future (past) directed geodesic 
normal. Let $ be the maximum of these distances. Let p 
be a point on a future (past) directed geodesic normal av 


a distance greater than È . Consider the function の with 


Z n 
っ ~ 4- - - + 1 E) d A + dy we 
respect to p over the compact surface li”. Let be tne 


Zz 
geodesic from p normal to H^ at the point q, wnereO has its 


maximum. ‘here must be a point conjugate to Hi? along A in 
But if there is no point conjugate to q along in ap, then 
A cannot be maximal by the first lemna. If however there 


is a point conjugate to q along À in qp, then tnere must be 


a longer Geodesic from q to p by tne second lemma. thus 


AN 


is not the geodesic of maximum length from H^ to p. ‘his i 
a contradiction which shows that the original assumption ta 
the space was non-singular must be false. 

his proof could also be used to show the occurrence 
of a singularity in a model with a non-compact Cauchy surf 
provided tnat the expansion of its normals was bounded away 
from zero and provided that the intersection of the Cauchy 


surface with all the time-like and null lines from & point 


was compact. 
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